A FAST ELEMENTARY ALGORITHM FOR COMPUTING THE 
DETERMINANT OF TOEPLITZ MATRICES 
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Abstract. In recent years, a number of fast algorithms for computing the determinant of a 
Toephtz matrix were developed. The fastest algorithm we know so far is of order log n + k^, 
where n is the number of rows of the Toeplitz matrix and k is the bandwidth size. This is possible 
because such a determinant can be expressed as the determinant of certain parts of n-th power of 
a related k X k companion matrix. In this paper, we give a new elementary proof of this fact, and 
provide various examples. We give symbolic formulas for the determinants of Toeplitz matrices in 
terms of the eigenvalues of the corresponding companion matrices when k is small. 



1. Introduction 



In this paper, we consider a.n n x n Toeplitz band matrix T„ with r and s super- 
diagonals as shown below: 



a-s+i 



T 



ao 



(^s+r 



\ 



ai 



V 



a^+r ■ ■ ■ as+1 ao J 



where r + s = k, as ^ 0, Ug+r 0, ai E K for i = 0, k and is a field. Without 
loss of generality we assume that s < k, as the determinant remains the same under 
taking the transpose, and our main concern is the determinant of r„. 

Finding fast algorithms to compute det(r„) is of interest for various applications. A 
number of fast algorithms computing det(T„) are developed recently (for tri-diagonal 
and pentadiagonal cases, see [3], [1], [5], [7j, [10] )■ Whenever r = s = 2, the author 
[3] gave an elementary algorithm computing det(T„) in 82y/n + O(logn) operations. 
In this paper, we both improved and generalized the algorithm given in [3]. Namely, 
the algorithm we give here works for any r > 1 and s > 1, i.e., for any k > 2. 
Moreover, it takes 0(|/c^ logg f + s^) operations to compute det(T„). The key part in 
this improvement is that the computation of det(T„) can be related to the powers of 
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the following kxk companion matrix C associated to T„: 



'1.11 



C 



-ai 
—do 

Is 

-Q3+1 



\ 



'(fc-l)x(fe-l) 



\ ■■• / 

\ as / kxk 

where /(fc-i)x(fc-i) is the identity matrix of size (A; — 1) x (A; — 1). The characteristic 
polynomial chc{x) of C is given by 

det(a;/ - C) = + ^x^-' + ■■■+ '^x^-^ + '^x^-^ + ■■■ + + " 



Note that the algorithm we give in this paper is the fastest known algorithm to 
compute det(T„), and was first given in [1]. Our method is elementary and more 
clear. More precisely, we give a new proof of the following theorem: 

Theorem 1.1 (Theorem 12.71 in [|2]). Let M be the upper left s x s suhmatrix of . 
For every integer n > k, det(T„) = (— l)"'*a" ■ det(M). 

After this theorem, our attention focuses on the computations of powers of C. Thus, 
we give brief discussion of various methods of computing C" in ^ 

Closed form formulas for det(r„) can be given in terms of the roots of chc{x) as 
previously described in [13]. In ^ we illustrate how this is possible for pentadiagonal 
Toplitz matrices, and give explicit formulas. 

The method given here is effective as long as k, the number of nonzero diagonals 
of T, is not close to n which is the number of rows of T. 



2. A FAST ALGORITHM FOR COMPUTING det(r„) 

In this section, we give an elementary algorithm for computing det(T„), which is 
the fastest known algorithm so far. 

First, we describe the outline of the algorithm as follows. We move the first s column 
vectors of T„ and make them the last vectors, successively. This gives a matrix P. 
If the column vectors of T„ are {Ci, C2, ■ ■ ■ , C„}, then the column vectors of P are 
{Cs+i, Cs+2, . . . , C„, Ci, . . . , CJ. Note that det(r„) = (-l)("-i)^ det(P). Then we 
multiply the first column by suitable terms and add to the last s columns of P so 
that the only nonzero entry in the first row will be a^. If the resulting matrix is P', 
det(P) is nothing but times the determinant of the cofactor P{ ^ of P'. We note 
that PI I is a matrix of size (n — 1) x (n — 1), and it is of similar form as P. Following 
the same procedure applied to P, we relate det(P{ J to the determinant of a matrix 
of size [n — 2) x [n — 2). Continuing in this way, the problem of computing det(Tn) 
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can be reduced to the computation of the determinant of a A; x /c matrix, and this 
matrix can be computed easily as it is n-th power of a kx k companion matrix. Next, 
we describe this algorithm in detail. To be precise, we introduce some notation and 
deduce some results. 

Let / : M.^ — y M.^ be a linear transformation given by /([xi, X2, ■ ■ ■ , XkY) = 

[X2 — Xi-^-^, X^ — Xi-T^, . . . , Xs+1 ~ XlTT} Xs+2 ~ 2^1-^, . . . , Xk — Xi ^ , — Xi^]*, 

where is the transpose of a vector v. Let F be a map sending a k x s matrix A to 
another kxs matrix F{A) by applying / to every column of A. That is, if the columns 
of A are {Ci, C2, . . . , CJ, then the columns of F{A) are {/(Ci), /(C2), . . . , /(C,)}. 

We define a sequence of k x s matrices (A)i>o recursively by setting A+i = -^(^0 
for every integer i > 0, and by taking the following initial value: 



/ Oq CLl ■ ■ ■ 0,3-2 \ 
O^s+l Oq tti ■ ■ ■ as-2 



s+r 



\ 



ao 

O-s+r / 



kxs 



Let {vi, V2, . . . , Vk} be the standard basis, consisting of column vectors, of M*^. 
Note that f{vi) = . . . , ^, ^, . . . , and /K) = v..^ for each 

i = 2, 3, . . . , k. Thus, the matrix C given in Equation (11. ip is nothing but the matrix 
representation of the linear transformation / : M.^ — > M'^. 



Lemma 2.1. Let the transformation F and the matrices Ai and C he as defined 
before. We can express F in terms of C as follows: F[Ai) = C ■ Ai. Therefore, for 
any integer i > 1 we have Ai = C ■ Aq. 



Proof. Note that C ■ A^ is nothing but Aj+i. Thus, the proof of the first part follows 
from the definition of F. Applying the first part successively to Aq gives F^{Aq) = 
C* ■ Aq. On the other hand, F'^{Aq) = Ai. This completes the proof. □ 
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Let 0(„_fc)xs be the zero matrix of size [n — k) x s. We set 

/as \ 

o-s-i as 



Br, 



ao 

Cls+r ag^r—i 



\ 



/ 



p . 



Bn 



A, \ 



{n-k)xs f ^ 



nx(?i— s) 

Next, we relate the determinants of T„ and P„o- 



Lemma 2.2. For any n > k and s>0, det(r„) = (-l)("-i)Met(F„,o)- 

Proof. We move the first s column vectors of T„ and make them the last vectors, 
successively. This gives the matrix P. More precisely, if the column vectors of T„ are 
{Ci, C2, ■ ■ ■ , Cn}, then the column vectors of P are {Cs+i, Cs+2, ■ ■ ■ , Cn, Ci, ■ ■ ■ , Cs}- 
Note that the matrix P is nothing but Pnfl. Therefore, det(T„) = (— det(Fn,o)- 

□ 

The following theorem is a generalization of [3l Theorem 2.3]: 

Theorem 2.3. Let s and k be as before. For every integers n > k and i such that 
< i < n — k, we have 

det(r„) = (-l)("-i)X ■ det(P„_,,). 



In particular, 



det(r„ 



-1) 



{n-~l)s n—k 



det(Pfc, 



n~k I 



Proof ByLemma[221 det(r„) = (-l)("-i)Met(F„,o)- Therefore, we are done ii n = k 
or i = 0. For the rest of the proof, we assume n > k and z > 0. 

For any integer i with 1 < i < n — k, we multiply the first column of Pn-i+i,i~i 
by and add it to its (n — i — s+j + 2)-th column for each j = 0, 1, . . . , s — 1, 
where bj is the + l)-th entry of Let Ri be the resulting matrix. Clearly, 

det(P„_i+i,i_i) = det(i?i). The only nonzero entry in the first row of Ri is the (1, l)-th 
entry with value a^. Note that (1, l)-th minor of Ri is nothing but Pn~i,i- Expanding 
the determinant of Ri using the first row gives det(i?i) = ■ det(Fn-i,i)- Thus, 
det(P„_j+i,i_i) = Os ■ det(F„_i,i). 

This gives that det(P„_o) = a\ ■ det(P„_i,i) for each integer i with < z < ri — fc. 
Then the result follows. 
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□ 

Theorem 2.4. For every integer n > k, we have 

det(r„) = (-l)("-i)^ar'= ■ deti[Bk,C^-' ■ Ao]). 

Proof. By the definition Pk,n-k = [Bk \ An-k]- On the other hand, An-k = C""'^ ■ 
by Lemma 12.11 with i = n — k. Then the result foUows from the second part of 
Theorem O □ 

We can improve Theorem 12.41 as follows: 

Theorem 2.5. Let M be the upper s x s suhmatrix of C"'~'^ ■ Aq. For every integer 
n > k, we have 

det(T„) = (-l)("-i)^ar' ■ det(M). 

Proof. The proof follows by similar arguments given in the proof of Theorem 12. 3[ but 
we need to introduce new notation for precise and shorter description. 
We have 

(2.1) det(T„) = i-lY^-'>a:-' ■ det([i?fe | 

by the second part of Theorem 12.31 and the proof of Theorem 12.41 

Let Ui denote the lower right [k — i x r — i) submatrix of B^, and let Vi denote the 
upper k — i X s submatrix of An-k+i for i = 0, 2, . . . , r. Note that [Ur | K] = the 
upper s X s submatrix of An-s, and \Uq \ Vq] = [Bk I ^n-fc]- 

For any integer i with < i < r, we multiply the first column of \Ui \ V^] by ^ 
and add it to its {r — i + j + l)-th column for each j = 0, 1, . . . , s — 1, where hj is 
the (1, j + l)-th entry of Vi. Let Si be the resulting matrix. Clearly, det([f7j | V^]) = 
det(S'i). The only nonzero entry in the first row of Si is the (l,l)-th entry with 
value Og. Note that (1, l)-th minor of Si is nothing but [f/j+i | V^+i]. Expanding the 
determinant of Si using the first row gives det(5'i) = Os ■ det([[/i+i | V^+i]). Thus, 
det([/7i|y,]) = a,-det([[/,+i|V^,+i]). 

This gives that det([[/o | Vq]) = ■ det([L^r I K])- That is, det{\Bk \ An-k\) = al ■ 
det(K)- Combining this with Equation (El]), we see that det(T„) = (-l)^"-^)*^-" • 
det(M), where M is the upper s x s submatrix of An-s- On the other hand, An-s = 
(jn-s . Lemma 12.11 Then the result follows. □ 

Note that Theorem 12.51 outlines a fast algorithm to compute det(T„). However, we 
look for even simpler formula, which is possible if we find an answer to the following 
question: 

Can we get a simpler expression for the determinant of the upper s x s submatrix 
of C^-' ■ AqI 

Lemma 12.61 below provides a fairly simple expression for this determinant. First, 
we give some well-known facts that we will use. 
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For each integer i = 1, 2, ... A;, suppose Mj „ satisfies the recurrence relation 

_ ao ds+l Cls+r-l Cls+r 

•^n+k — -^n+k—l ' ' ' -^n+k—s -^n+r—l ' ' ' -^n+l -^n; 

and that their initial values are given by the first equality below, where /fexfc is the 
identity matrix of size k x k. Then the powers of C are given by the second equality 
below: 

/ Ul,k-1 ■ ■ ■ Uifi \ ( Xii^„+fc_2 • • • Ml,n \ 



-^fcxfe, C*" 



\ Uk^k-\ Uk,k-2 • ■ ■ Ukfi ) 



U2,n+k-l U2,n+k-2 ' ' ' ^2,1 



\ Uk,n+k-l Uk,n+k-2 ' ' • Uk,n / 



Lemma 2.6. For every positive integers n and s with n> s the determinant of 
the upper s x s submatrix of C"'~^ ■ Aq is {—asY times the determinant of the upper 
left s X s submatrix of C"'. 

Proof. Let ci^m be the (/, m)-th entry of the matrix C"""* • Aq of size k x s. Since ci^m 
is Z-th row of C"~* times m-th column of Aq, we have 

m r 
j=l i=l 

Using the recursive formula of ui^n+k-m, we can express the above equality as follows: 

s—m 

Now, we note that m-th column of C"~* • Aq contains hnear combinations of its last 
s — m columns, and that the determinant does not change by adding a multiple of a 
column to some other column. 

Let Mq be the upper s x s submatrix of C™~^ ■ Aq. For every t = 0, ... s — 2, 
suppose Mf+i is the matrix obtained from Mf by multiplying its (s — t)-th column by 
— and then adding to its {s — t — j)-th. column for each j = 1, 2,..., s — t — 1. 
Note that det(Mo) = det(Ms_.i) and that (/,m)-th entry of Ms-i is —asUi^n+k-m- 
That is, det(Ms_i) is {—dsY times the determinant of a matrix with (/,m)-th entry 
ui^n+k-m, which is nothing but the upper left s x s submatrix of C". This completes 
the proof. □ 

Next, we give the main result of this paper. Namely, the computation of det(T„) 
is basically reduced to the computation of n-th power of the k x k matrix C, and the 
computation of the determinant of an s x s matrix: 

Theorem 2.7. Let M be the upper left s x s submatrix of C"'. For every integer 
n>k, 

det(Tn) = (-l)"X-det(M). 
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Proof. The result follows from Theorem 12.51 and Lemma 12.61 □ 

Since M is s x s matrix, its determinant can be computed in 0{s^) operations. 
In fact, the power 3 here can be taken less. More costly part is the computation 
of C", which can be done in 0(|fc^log2 f ) operations as explained at the end of ^ 
Therefore, det(T„) can be computed in 0(|A;^log2 | + s^) operations. 

Note that Theorem 12.71 is very similar to Prop 2.1]. However, the proof we gave 
here is more elementary and clear. 

We can use the algorithm described in this section to compute the characteristic 
polynomial of T„ if we start with T„ — XI rather than T„. Therefore, Theorem 12. 71 can 
be extended as follows: 

Theorem 2.8. Let C\ be the matrix obtained from C by replacing ao by oq — A. For 

every integer n > k, we have 

det(T„ - XI) = ■ det(MA), 

where I is the n x n identity matrix and M\ is the upper left s x s submatrix of C^. 

3. Computing 

In this section, we briefly describe the several ways to compute C" from C. 
We set C° = /, where / is the identity matrix of size k x k. Then we use = 
(jn/2Qn/2 j£ ^ -g gyg^^ g^^^j ^gg (jn _ (j(jn-i ^ ^^lis way, C" cau bc 

computed in O(logn) times the number of operations to compute the product of two 
k X k matrices. If k is large, one should consider other methods of computing C". 

As shown in the previous section, can be expressed in terms of k homogeneous 
linear recurrences Uij, where i G {1, 2, ■ ■ ■ , k}. If we compute the roots of the 
characteristic polynomial of C, chc{x), we can express each of these recurrences in 
terms of these roots since they satisfy the recurrence relation given by chc{x) and 
their initial values are known (see [HI section 7.2.9] for solutions of such recurrences). 
We should note that the computation of roots of a polynomial of degree higher than 
five can be a difficult task. 

Another method is to use Jordan canonical form of C. We explain the advantageous 
steps in this approach. First, we note that the characteristic polynomial of C, chc{x), 
is the minimal polynomial of C (see [TTl page 325]). Over C, chc{x) factors into 
linear terms. Then C have a Jordan canonical form J such that C = V~^JV, for 
the Vandermonde matrix V as shown in Equation (13. ip (if any of the eigenvalues has 
multiplicity more than one, then V will be the corresponding confluent Vandermonde 
matrix). Then C" = V^^J^V. This along with Theorem 12.71 implies that det(r„) = 
(— l)"'*a" ■ det(M), where M is the upper left s x s submatrix of V~'^J"'V. Moreover, 
we have the following useful information about the part J": 

If J = diag(Ji, J2, . . . , Js) with s < k and Jj is Jordan block matrix for I < i < k, 
then J" = diag( J", J^, . . . , J") for every integer n > 1. 
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Let J(i)(A) be a Jordan block matrix, part of a Jordan matrix J, corresponding to 
an eigenvalue A. As shown in Equation (I3.ip . J(i){X) is an upper triangular square 
matrix of size i x i and having A's on its diagonal, I's on the next diagonal and O's 
elsewhere. 



(3.1) J(,)(A) 



/A 1 
A 1 

A 



1 

A/ 



V 



( A^i \\ 



, fc-i 



V A 



fc-i 

k 



k-2 



k-2 



k-2 



Ai 1 \ 
A2 1 

Afe 1 / 



kxk 



Then one can show that [12, Exercise 7.3.27] its n-th power is as follows: 



^(•)(A) 



\ 



A" 



nA"-i 
A'^ 



n{n-l) \ n-2 

nA"-i 
A" 



where J(")(A) has n + 1 nonzero diagonals (or less if the dimension of J(i)(A) is less), 
and the j-th diagonal from the main diagonal has entries (")A"~-^. 

Note that various methods of matrix exponentiation are explained in the article 
[8]. One can use the modified versions of those methods to compute C". The reader 
can find more information on the critique of those methods in the article [8] and the 
references therein. 

One can use combinatorial arguments to derive combinatorial formulas for the en- 
tries of any power of C (see [2] and [6]). Although such combinatorial formulas are 
very interesting, they don't lead to fast algorithms to compute powers of C. 

Assuming that the roots of chc{x) are known, one can compute in ©(/clogg^) 
operations by using the method proposed by W. Trench [13J. D. Bini and V. Pan 
[H Page 436] improved Trench's result. Namely, without assuming that the roots 
of chc{x) are known, one can compute in (^(Ifc^loggf) operations. If the field 
that C defined over supports Fast Fourier Transform, Bini and Pan showed that C" 
can be computed in 0(A;log A;log2 |) operations (see [U Page 436] and the references 
therein) . 
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4. Determinants of tridiagonal Toeplitz matrices 

In this section, we consider the tridiagonal matrices, i.e., we have s = 1, r = 1, and 
so A; = 2, 

(ah \ 

c a 



T 

■J- 71 



■•. h 

c a 



and C 



■f 1 
-% 



/ Jixn 



where 6 7^ 0, c 7^ 0. In this case, the characteristic polynomial of C is det(a;/ — C) 



We note that 



where both Un and Wn satisfy the recursive 



relation Xn+2 = -Y^n+i — for each n > 0, and they have the initial values uq = 0, 
ui = 1, Wo = I, wi = 0. We recognize that is nothing but the Lucas sequence 
of the first type. That is, u„ = Un{P,Q) for each integer n > 0, where P = ^ and 
Q = |. Because f/„+2(-P, Q) = PUn+i — QUn for each integer n > 0, Uo{P, Q) = and 
Ur{P,Q) = l. 

Applying Theorem 12.71 gives 

det(T„) = (-ir6'^^„,+i(^,^). 

Hence, existing formulas for Un applies to det(r„). For example, using the binet 
formula for Un+i we obtain 



det(T„ 



- 46c 



and using the recursive formula for Un+i we obtain 

det(r„+2) = a ■ det(T„+i) — be ■ det(T„), for n> 2, 
where det(Ti) = a and det(T2) = — be. 

5. Determinants of pentadiagonal Toeplitz matrices 

In this section, we give explicit formulas for determinants of pentadiagonal matrices. 
In this case, we have s = 2, r = 2, and so A; = 4, 

/ a b e \ 

/ 1 \ 

-2010 



d a 
e ' • . 







a b 
e d a 



and C 



-f 1 
V -f / 
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where c 7^ 0, e 7^ 0. The characteristic polynomial of C is chc{x) = det(a;/ — C) = 
+ -x^ + -x^ + -x + -. 

c c c c 

For each integer n > 1, C" is expressed in terms of Un, Vn, tn and Wn, each of which 
satisfies the recursive relation Xn+A = 



lXn+3 - ^Xn+2 - -Xn+1 - 'Xn, and that their 



initial values are given by the second equality below: 



/ U„+3 Un+2 



t 



n+3 



t 



n+2 



Vn+l 
tn+l 



Un \ 

Vn 



, and 



\ Wn+3 Wn+2 Wn+l Wn J 



U2 
V2 
t2 



Ui 



to 



\ W3 W2 Wi Wo J 



/ 1 \ 

10 

10 

1 



Applying Theorem O gives det(T„) = c" ■ {un+3Vn+2 - Un+2Vn+3)- 

As done in the previous section, one can use the formulas for the homogeneous 

linear recurrence relation Un and Vn to derive a formula of det(T„) in terms of the 

roots of the characteristic polynomial of Un and Vn, which is nothing but chc{x). 

Instead, we used the Jordan canonical forms in our computations below: 

Note that C is an irreducible Hessenberg matrix, and that both the minimal and 

the characteristic polynomials of such matrices are the same (see pT| page 325]). 

Alternatively, the reader may use a computer algebra system such as Mathematica 

[9]) to check that the minimal polynomial of C can not be of degree 1, 2 or 3. This 

makes some restrictions on the types of Jordan canonical forms, since the multiplicity 

of an eigenvalue A in the minimal polynomial is the size of the largest Jordan block 

corresponding to A. We have the following five cases. We used Mathematica [9] for 

the details of these computations. 

Case I: chc{x) = {x — Ai)(x — A2)(a; — A3)(x — A4), where the eigenvalues Aj's for 

i = l,---4 are distinct. Then, J = diag( J(i)(Ai), J(i)(A2), J(i)(A3), J(i)(A4)), and we 
E 

have det(T„) = c"— , where 

E={X2- A3) (Ai - A4) ((A2A3)"+' + (AiA4)"+') - (Ai - A3) (A2 - A4) • 
(5.1) ( (AiA3)"+' + (A2A4)"+' ) + (Ai - A2) (A3 - A4) ((AiA2)"+' + (A3A4)"+') , 
D = (Ai - A2) (Ai - A3) (Ai - A4) (A2 - A3) (A2 - A4) (A3 - A4) . 



Case II: chc{x) = {x — Ai)^(x — X2){x — A3), where the eigenvalues Ai, A2 and A3 

are distinct. Then, J = diag( J(2)(Ai), J(i)(A2), J(i)(A3)), and we have det(r„) 
where 



(5.2) 



E = Aj+"(At+"(A2 - A3) + Af "(Ai(-(2 + n)X, + A2 + nAa) + ((3 + n)Ai 
- (2 + n)A2)A3)) + A^+"(A^+"(A2 - A3) + A}+"((2 + n)A? 
+ (2 + n)A2A3 - Ai((3 + n)A2 + A3 + nX^))), 

D=(Ai-A2)'(Ai-A3)'(A2-A3). 
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Note that formula of ^ using Equation (15 ■2p is the hmiting value of ^ using Equa- 
tion (15. ip as A4 — )■ Ai- 

Case III: chc{x) = (x — Ai)^(a; — A2)^, where the eigenvalues Ai and A2 are distinct. 

Then, J = diag( J(2)(Ai), J(2)(A2)), and we have det(T„) = c" — , where 

E = Af 2" - (2 + n)2(AiA2)^+"(A? + A^) + 2(3 + 4n + n2)(AiA2)'+" + A^+^n^ 

(5.3) 

D = (Al-A2)^ 

In this case, the formula of ^ using Equation (I5.3P is the limiting value of ^ using 

Equation (15. 2p as A3 — > A2. 

Case IV: chc{x) = (x — Ai)'^(x — A2), where the eigenvalues Ai and A2 are distinct. 

E 

Then, J = diag( J(3)(Ai), J(i)(A2)), and we have det{Tn) = c"-^, where 



E = i2 + n)\^(il + n)(A?+" - A3+") - (3 + n)A?+"A2 + (3 + n)AiA2+- , 
(5.4) 

= 2 (Ai - A2)' . 

Note that the formula of ^ depending on Equation (15. 4p is the limiting value of ^ 
relying on Equation (15. 2p as A3 — )■ Ai. 

Case V: chc{x) = (x — Ai)^. Then, J = J(4)(Ai), and we have 

det(r„) = _(n + 3)(n + 2)\n + 1)A^. 

Note that this formula of det(r„) is the limiting value of det(T„) given in case IV as 
A2 — > Ai. 

Next, we provide examples for each of these five cases. 

Example I: If a = 101, b = -17, c = 1, d = -247, e = 210, then we have 
chc{x) = {x — 2){x — 3){x — 5){x — 7), and for any integer n > 4, 

det(T„) = Y^(-6 • 10"+^ + 5 • 15"+2 + 6"+^ - 6 ■ 21"+^ + 5 ■ 14"+^ + 35"+2). 

Example II: If a = 17, 6 = —7, c = 1, d = —17, e = 6, then we have chc{x) = 
{x — l)^(x — 2){x — 3), and for any integer n > 4, 

det(T„) = ^(2"+2(2n + 3) - 3"+2(2n + 5) + 6"+' + 1). 

Example III: If a = 37, 6 = —10, c = 1, d = —60, e = 36, then we have 
chc{x) = (x — 2)^(x — 3)^, and for any integer n > 4, 

det(r„) = 4"+2 + 9"+' - {n{n + 4) + 16)6"+^ 

Example IV: If a = 30, 6 = —9, c = 1, d = —44, e = 24, then we have chc{x) = 
{x — 2Y'[x — 3), and for any integer n > 4, 

det(T„) = 2"-i(n + 2){T+^{n - 3) + 2"+^{n + 7)). 
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Example V: If a = 24, b = —8, c = 1, d = —32, e = 16, then we have chc{x) = 
(x — 2)^, and for any integer n > 4, 



It is straightforward to implement the algorithm we outlined to obtain similar closed 
form formulas for other small values of r, s (and so k) by using Mathematica [S] or 
some other computer programs having symbolic capabilities. However, the formulas 
for A; > 6 are not short enough to include here. 
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